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>- 1 Abstract 

o 

I In this study the determinant of the average quadratic error matrix is used 

as the measure of state estimation efficiency. Minimizing this quantity gives 
us the optimal measurements in different scenarios. We present applications 
when von Neumann measurements or a single POVM are used, when there is 
no known information or a part of the parameters of the state is given. Under 
some restrictions the optimality is found for n-level systems. The optimal mea- 
surements have some complementary relation to each other and to the available 
data, moreover, symmetric informationally complete systems appear, containing 
a new, conditional version. 
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1 Introduction 



State estimation or tomography is a fundamental problem in the field of quantum in- 
formation theory and it can be considered as one of the foundational issues of quantum 
mechanics [HE]- The topic of quantum tomography consists of methods which recon- 
struct the state of the system under investigation using repeated measurements of a 
set of observables. The problem may be traced back to the seventies [11], the interest 
in a thorough mathematical analysis of the quantum state estimation procedures has 
been flourishing recently [21 [TUl H21 125] . 

In statistics the accuracy of the estimation can be quantified by the quadratic error 
matrix. The matrices are typically not comparable by the positive semi-definiteness, 
hence if different estimation schemes are compared, the determinant of mean quadratic 
error matrix can be used instead. This approach can be found in references [T9| 120] , 
their main result was that the complementary von Neumann measurements are optimal. 
A more general context appears in [3], which is much closer to our approach. The 
concept of complementarity was extended to quantum subsystems in [2TJ and the case of 
two qubits has a detailed analysis in [IJ1I2I]. When all parameters of the density matrix 
are obtained from a single measurement, then a symmetric informationally complete 
POVM appears [26J. A similar result was obtained earlier by Wootters and Fields (29] 
for von Neumann measurements and by Scott [28J for POVMs, but optimality had a 
different formulation in both cases. 

Some a priori information about the state can be given in various ways, for example 
one can minimize the relative entropy to an a priori given state, resulting in a bias 
toward that state [15]. The most popular subject in this field is state discrimination: 
when we know that the system is one of several given states and we should figure out 
which one [5]. Beside knowing the possible states we can have an a priori probability 
distribution on the true state, too. This idea was used in [B] to obtain the optimal 
phase estimation. The given states do not construct a discrete set, instead, they are 
searching among all pure states. Another continuous subset is used in [12], where we 
know that the state to estimate is a convex combination of two given states. In our 
setup we know that the state is on a given subset of the whole state space (this can 
be the whole state space itself) and assume an arbitrary unitarily invariant a priori 
probability distribution to be known. 

We use the determinant of the average quadratic error matrix as a figure of merit 
in state estimation setup; the average is taken over the unitarily invariant states. This 
quantity is relies on purely classical statistics, therefore it is much easier to compute 
than, e.g., quantum Fisher information, but still reproduces the well-known quantum 
phenomena and even more. So it gives us a useful tool to find the optimal measurement 
setup for different quantum tomography problems. We examine the assumption that 
the unknown state is known on some complementary subalgebras, which practically 
means that one should restrict the are of integration. For our calculations we use a 
parametrization for both the density operators and the observables (we expand them 
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in a fixed basis that is orthonormal with respect to the Hilbert-Schmidt inner product), 
but we can see that the results are independent of this parametrization. To gain ana- 
lytic results, some sort of symmetry is recommended but by using numerical methods 
we could overcome that and solve difficult problems efficiently. 

Section 2 is a very short introduction, quantum states, measurements, state esti- 
mation, complementarity and symmetric measurements are overviewed. In Section 3, 
we give the motivation for using the determinant of the average covariance matrix by 
guiding through a simple example, and then give some useful remarks about this quan- 
tity. Then the optimization problem is solved for arbitrary dimensions under some 
restrictions. In Section 4 the complementarity of optimal von Neumann measurements 
is obtained. In Section 5 a single POVM is used for the measurement, then a symmet- 
ric informationally complete POVM appears. Finally, we investigate the asymmetrical 
case for POVMs too (Sec. 6) and draw conclusions (Sec. 7). 



2 Basic concepts 



In this section, we give a brief summary of the basic concepts appearing in the present 
paper (for further details, see [HE]). 

p G M n (C) is an n-dimensional quantum state if 



p > and Tr p = 1 . 



(1) 



We will use a basis of self-adjoint matrices to obtain a parametrization of these 
density matrices. In the qubit case, we use the Pauli matrices as a basis, i.e., the Bloch 
parametrization: 



1 



(I + 9 -a) 



i=i 



1 + # 3 01 - 1#2 

e l + i6 2 i - # 3 



(2) 



In this case the whole state space of quantum states can be described using the 
Bloch vector: 9 = (0i,9 2 ,0 3 ) and the conditions in flTJ are simply converted into 
01 



91 



el < i. 



For the n-dimensional case we use the following parametrization: 

1 



P 



n 



-(1 + 6 -a) 



where 6 G 



is 



n 



the generalized Bloch vector, 6 ■ a = Y2j QjCj an d {cj : 1 < j < 
1} are generalized Pauli matrices, which form an orthonormal basis on the self- 



adjoint traceless matrices: Oi 
3-dimensional case are the multiples of Gell-Mann matrices. 



<rt, Tra, 



and Tr<7j<7j 



5i t j. An example for the 
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A positive operator-valued measure (POVM) can in this framework be described as 
a set of matrices: 

k 

E = {E 1 ,E 2 ,...E k } with ^Ei = I,Ei>Q. 

i=l 

Von Neumann measurements can be described as a POVM with two components and 
projective elements: {Pi, / — Pi}. 

Having a quantum state p and a POVM E, the probability of obtaining an outcome 
related to the POVM element Ei is: 

Pi = Ti{pE i ). (3) 

If we have m samples of the state p and we measure P on each copy independently, 
then let z/j be the relative frequency of measuring the ith outcome. If m is a large 
number, then z/j will be close to p iy and if we know the POVM (because we set it), 
then we can make an estimation on p from ([3]) by using z/j instead of pi. A simple state 
estimation problem is described in detail in Section [3j 

The heuristic concept of complementarity was born together with quantum the- 
ory. A mathematical definition is due to Accardi [1J and Kraus [3]. Let % be an n- 
dimensional Hilbert space. Let the observables A and B have eigenvectors ei, e 2 , . . . , e„ 
and fx, f 2 , ■ ■ ■ , fn which are orthonormal bases, they are complementary if 

{{eiJ^^l (l<*,J<n). (4) 

If this condition holds then the two bases are also called mutually unbiased. An 
overview about complementarity can be found in [22J and details are in [2TJ 124] . 

Complementarity can be generalized to the case of POVMs. The POVMs {Pi,P2, 
. . . , Pfc} and {Pi, F 2 , . . . , F m } are complementary if 

Tr EiFj = -Tr Ei Tr F j (l<i<k, 1 < j < m). (5) 

This is equivalent to the orthogonality of the traceless parts: 

TrPi TrP 
Ei -I JL F 3 3 -I, 

n n 

so we will use the expression quasi-orthogonal, when we use this property of comple- 
mentary operators. 

The concept can be extended to subalgebras A\,A 2 C M n (C). Ai and A 2 are 
complementary if A\ 6 A\ and A 2 G A 2 are quasi-orthogonal (that is, the traceless 
part of matrices are orthogonal). 

An orthonormal basis corresponds to the subalgebra of operators which are diagonal 
in this basis. Such a subalgebra is maximal Abelian, it will be called M-subalgebra. 
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A symmetric informationally complete POVM {Ei : 1 < i < k} of an n-level system 
is described by a set of projections Pj (1 < i < k) such that 



where k = n 2 , A = n, fi = l/(n + 1) and all projections are rank-one, see [26] [2?j. The 
existence of a symmetric informationally complete POVM (SIC POVM) is not known 
for every dimension n. 

In this paper k < n 2 case will also appear, when some parameters are already known, 
we need to estimate fewer parameters, i.e., fewer elements in POVM are needed to make 
a full reconstruction. If they fulfill condition ([6]) with some constants and the P r s are 
complementary to the known subalgebra, then we call them conditional symmetric 
informationally complete POVM toward the known n 2 — k parameters, because they 
are informationally complete only if we use the condition that we want to discriminate 
the states only in the unknown directions. They include the case of SIC-POVMs as a 
special case when no parameters are known. The existence of such a POVM can be a 
fundamental question in different quantum tomography problems. 

3 The determinant of the average covariance ma- 
trix 

Suppose that we have a qubit and some of the three parameters 8\, 82 and 63 are known 
and the unknown parameters will be estimated. For example, let us assume that 8\ 
and 82 are known and we want to estimate #3. The assumption means that the reduced 
state is known on the M-subalgebras generated by U\ and 02, respectively. They are 
complementary subalgebras and the M-subalgebra generated by er 3 is complementary 
to both. A projection of a von Neumann measurement 



is used to estimate #3 from the result of several measurements, where Af + A| + A§ = 1. 
The expected value is 



k 




(6) 



i=l 



E = \{I + \-a) 




l + ^Tra.P 



3 



) 



^(i + (M» 



and we have 
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Let us denote by v the random outcome of the measurement of E, the expected value 
is p = Tr pE and it is Bernoulli distributed. The natural unbiased estimate of (9 3 is 



The variance is 




(7) 



Var (9 3 ) = -^H^ = . (8) 



We want to argue that the optimal estimate corresponds to the projection Ai = 
A2 = 0, A3 = ±1, i.e. measuring in the unknown direction. This estimator has the 
variance 1 — 6\. However, the inequality 

1 - (A,0> 2 n o 

A|^- 3 (9) 

is not true in general, for example A = 9 is possible for pure states and then the left- 
hand side is zero. So if we want to minimize Var (#3), it is not true that ; ^p- is the 
best observable. 

Since 9 1: 9 2 are given, the only unitarily invariant vector to (#1, 9 2 , 9 3 ) is (9i, 9 2 , —#3), 
if we take the average of the left hand side of Q on these states we obtain 

1 1 — (Ai#i + X 2 9 2 + A 3 # 3 ) 2 1 1 — (X\9 1 + X 2 9 2 — A 3 6> 3 

Var r 



2 A3 2 A3 



We can solve the minimization problem 

Var (0 3 ) ) -> min, Xj + A 2 + A 2 = 1 



using the Schwarz inequality, and get that in the optimal case Ai = X 2 = 0, A3 = ±1. 
Thus, the measurement of -^rp- is optimal indeed, yet only in the average sense. 

Note that in (J7|) v has Bernoulli distribution, with p(l — p) variance. If m mea- 
surements are performed, we will have v = (p\ + v 2 + . . . + u m )/m, where is the 
outcome of the z-th measurement. Then the variance will be mp ^~ p ^ = —p(\ — p), so 
we can conclude that if we have multiple measurements the variance changes only by 
a constant factor, hence the minimization problem remains unaffected. 

In the following sections, we will always use the same kind of naturally derived 
estimators; one can prove that those are unbiased and efficient. We should also always 
bear in mind that if we have only a few measurements, we can easily get a Bloch- 
vector as a result of the estimation, which does not satisfy the positivity condition 
in ([T]). But in practical cases, we have many measurements, hence the law of large 
numbers ensures us that the estimator will give a state near the original one, and the 
probability of getting a physically impossible state (a density matrix with at least one 
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negative eigenvalue) will converge to zero exponentially, due to the large deviation 
theory. That is why we do not have to deal with any physical restrictions, even if we 
calculate the variance only for one measurement in each case. 

If we have multiple parameters, we can still get a multidimensional estimator from 
the measurement outcomes, but in that case, a covariance matrix will take the role of 
the variance. The average covariance matrix can be calculated by taking the piecewise 
average of the elements, but since the matrices are not comparable in general, we 
need an M n (C) — > R function which gives us a number as the objective function to 
minimize. After examining some candidates we concluded that the best choice is to 
take the determinant of the matrix and that this needs to be done after the averaging 
process, so we have the following optimization problem: 

det (Gov (fl) ^ — > mm, 

where 9 is the estimator of the unknown parameters. This way we obtain an easily cal- 
culable quantity, and by solving the minimization problem, the optimal measurements. 

If we want to take the average of this covariance matrix on the whole state space, 
we have to give a unitarily invariant measure on the state space. On the other hand, 
if we integrate on the unitarily invariant states of an initial state with respect to 
the normalized Haar measure and the minimum is independent of the initial state, 
we get optimal measurements independently from the (unitarily invariant) probability 
measure. Thus, we can obtain a stronger result than we would by using a particular 
measure. 

Finally, let us note that using a POVM instead of von Neumann measurements is not 
a crucial difference. If there are n unknown parameters, we need to have n independent 
positive operators in the POVM. We only have the restriction = I, so n + I 

operators are sufficient and we can construct the estimator from the outcomes related 
to the first n components. The only difference is that in this case, the n outcomes 
would not be independent, off-diagonal elements appear in the covariance matrix, too. 



4 The optimal von Neumann measurement 

The n-level quantum system is described by the algebra M n (C). Consider the following 
decomposition 

M n (C) = CI®A®B, 

where A and B are linear subspaces and orthogonality is defined with respect to the 
Hilbert-Schmidt inner product (A, B) = Tr A*B. 

A state has the density matrix 

p = I/n + p A + p B . (10) 
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We assume that the component p B should be estimated. Let the dimension of B be 
k. The positive contractions E 1 ,. . . ,E k are used for independent measurements (on 
several identical copies of the n- level system): A measurement corresponds to the 
POVM {E l , I — E' 1 }. These operators have the expansion 

E i = ei I + E\ + E l B (l<i<k). (11) 

The expectations are 

p t := Tr pET = e t + Tr p A E A + Tr p B E B (1 < i < k). 

We fix an orthonormal basis Fi,...,F k in B. The unknown component has the 
expansion 

Pb = 0\Fi + . . . + 9kF k , 
where 6 — (9i, . . . , 6 k ) are the parameters to be estimated. Similarly, 

E B = e iX F x + e i2 F 2 + ... + e ik F k 

is an orthogonal expansion. 

The estimates §i are solutions of the equations 

k 

Ei = ei + Tr p A E\ + ^ e^Oj, 

3=1 

where £j is the random result of the ith measurement, 1 < i < k. In another form 







"ei" 




■Trp A E\- 




"en • 


• ei fc " 










+ 




+ 






, £ k. 








. Tr p A E A . 




.e fc i • 







or in a different notation 



1 n 



ei + Trp^E^ 



+ T9\ 



le k + Trp A E A ] 

where T is the k x k matrix from the previous formula. Therefore, 

9* = T- 1 ' 

We have 

§t_gt = T -i 



( 


"ei" 











( 


"ei" 




"Pi" 


) 













The quadratic error matrix is the expected value of 



i k 



[ei -Pi, ... ,e k - Pk \cr-^' - u - ■ ■ ' 

Ek -Pk. 

Due to the independence of the measurements the expected value of 

r 1 k 

(e i -p l )(e j -pj) 

L -l«,J=l 

will be diagonal. 

The random variables £,■ have Bernoulli distribution and their variance is 



We want to take the average: 



^1 — Tr pE i )TrpE i . 



(l-TrpE^TrpE^pip), 



where we integrate on the unitarily invariant states, and p is the corresponding nor- 
malized Haar measure. For the sake of simplicity, assume that the operators E l have 
the same spectrum. Then the integral is constant: it does not depend on the actual 
Ei, so the average of the quadratic error matrix is 



cT'^T-y. 

The determinant is minimal if the determinant of the matrix T is maximal. Geomet- 
rically, the determinant is the volume of the parallelepiped determined by the row 
vectors. To maximize the determinant, the row vectors should be long. This implies 
that E\ = in (TIT]) , since otherwise we could project E l on to CI © B, and still have 
the same elements in T. So we have 



E l = e + e ll F 1 + ... + e ik F k 

and then the determinant of the matrix 

1/n 
e en • 



(1 <i < k), 



e k i 





eifc 
e«c . 



is det T. The angle of the first row and any other row is fixed. To have a large 
determinant the rows of T should be orthogonal. In this case, the operators E 1 , . . . , E k 
are quasi-orthogonal. 
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Theorem 1 If the positive contractions E 1 , . . . , E k have the same spectrum, then the 
determinant of the average of the quadratic error matrix is minimal if the operators 
E 1 , . . . ,E k are complementary to each other and to A. 

Example 1 If two qubits are given and the reduced states of both qubits are known, 
then the ideal state estimation is connected to the observables 

011, (T22, 033, 012, 023, 031, 013, 021, 032, 

where — Oi® Oj. □ 



5 The optimal POVM 

Once again we are in the n-dimensional case, but in this section we will use POVMs 
for state estimation. Assume that A = {0}, which means that all parameters are 
unknown. Then 

E l = e l I + E l B = e l (I + fi-a) (1 < i < n 2 ), 

where £ R, /, G R n2_1 , /j • a = . fijOj and {aj : 1 < j < n 2 — 1} are generalized 
Pauli matrices, Trcr; = 0, Trof = 1, Tr<7j<x,- = 0. The positivity condition of is not 
known, but we have a necessary condition: 

Lemma 1 // the n x n matrix I + g ■ a is positive, then J2 k g\<r? — n. If J2 k g\ = 
n 2 — n, then I + g ■ a = nP with a projection P of rank 1. 

Proof. A = g ■ a is self-adjoint, Tr A = 0. Let Ai, A 2 , . . . , A n be the eigenvalues of A. 
Then 

n 2 - 1 n 

J29l = ^A 2 = J2^ 

k=l t=l 

Since Ylt=i — from Tr A = and A t > — 1 from I + A > 0, we have the upper 
bound. Namely, J^tLi ^ s maximal, if Ai, A2, . . . , A n is a permutation of the numbers 
— 1, — 1, . . . , — In this case 1 + A has eigenvalues 0, 0, . . . , 0, n, so it is a multiple 
of a projection. □ 

The computations are similar to the von Neumann case. The probabilities of different 
outcomes are 



Pi 




ei 




0i 








+ T 




-Pd. 




_e<i_ 




.o d . 



where d = n 2 — 1 and T — (ei, e 2 , . . . , ed) T (fi, /2, • • • , fd) dyadic matrix. 
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If Vi, u 2 , . . . I'd are the relative frequencies of the outcomes (of different measurements 
on identical copies), then the estimator is 



"01" 






- ei~ 




= T- 1 






Jd. 




yd 


- e d _ 



(12) 



The mean quadratic error matrix is 

V(6) = T _1 W (T -1 )* 

with W denoting the covariance matrix of the random variables v\, V2, ■ ■ ■ ^d- These 
are multinomially distributed, so 



W 



pi(l-pi) -P1P2 
-P1P2 ^2(1-^2) 



-piPd 



~P2Pd 



-PlPd 
-P2Pd 

Pd(l -Pd) 



To obtain the average mean quadratic error matrix, we integrate W on the rotation 
invariant states H with respect to the normalized Haar measure \x. We have to calculate 
two types of integrals: 

J -piPjd^O) = - J (ei + ei{fi, 0)) (ej + ejifj, 0)) dfj,(8) 

= - / eiejdfjiie) - eiej [ (f l ,9)d f i(9) - [ e i e j (f j ,9)d^(9)-e i e j [ (f t , 6)(f j , 9)dn,{9) 
Jh Jh Jh Jh 

= -titj + + - e i e j a(/ i , £). 

The integrals in the middle are zeros because of the symmetry: f H {v, 0) d/j,(9) = 0. 
The last integral comes from: 

^(f l ,9)(f J ,9)dfl(9) = jf (jlMt \ (YlUrnOn) = £ Qf ^M*)) /i^ 

The quantity J H 9jdfi(9) does not depend on £: it is a constant a depending on the 
domain of the integration; this way we obtained the stated formula. 

Similarly, 

/ Pi(l- Pi)d/i(9) = e^l - ti) - e*a(fi, /;), 
Jh 

so we can calculate the average quadratic error matrix: 



1-1 



e-'-l-aif,,^ -l- a (h,f 2 ) 
-l-a(hJ 2 ) e 2 1 -l-a(f 2 ,f2) 



-1 - a(fi,fd) 



-1 - a{f 2 ,fd) 



-1 - a{f u f d ) 
-1 - a(f 2 ,fd) 

e^ 1 - 1 -a(fd,fd 
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where F = (f u f 2 , . . . , f d ). 

The minimizer should be symmetric, hence = 1/n 2 , (fi, fi) = x and (fi, fj) = y, 
if i ^ j, z, j < n 2 — 1. Then the determinant has the form A/B, where 

A = (n 2 - a(x - y)) n ' 2 - 2 (l - a(x + (n 2 - 2)y)) 

and 

B = (x — yf- 2 ■ [x + [n 2 - 2)y) 

Therefore, we minimize 

A ( n 2 \« 2 -2 / j \ 

5 = l^" a J U+(- 2 -%" a J- (13) 

n 2 -l 

We can also calculate the length of f n 2 = — fi' 

i=i 

= - 1) • (x + (n 2 - 2)y) <n 2 -n (14) 

where the latter inequality is the condition for positivity, see Lemma [TJ On the other 
hand, we have the condition 

x < n 2 — n. (15) 

If both inequalities were sharp, then ffT3"|) would not be minimal, because we could 
increase both x and y with a sufficiently small e and then the value of (fT3l) would be 
smaller. If the equality holds in (I14p . then the second term of ( ITS"]) will be constant. 
So we want to have the difference of x and y as large as possible, hence x should be 
maximal (and then y is minimal). If the equality holds in fll5p . then from (I14j) we have 
y < —{n 2 — n)/(n 2 — 1), and on this domain f)13p has a negative derivative. So the 
minima are taken at 

x = n 2 — n, and y = 

Lemma [1] gives 

Ei=\(I + fi-<T) 

n z 

with some projections Pi and 

Tr = ^(1 + f l( r, I + f 3 a)-- 
We arrived at a system (J6j). The following statement is obtained. 

Theorem 2 If a symmetric informationally complete system exists, then the optimal 
POVM is described by its projections Pi as Ei = Pi/n (1 < i < n 2 ). 

For a qubit the existence of the symmetric informationally complete POVM is ob- 
vious, there are other known examples in lower dimensions, but in n dimensions the 
existence is not known. 



n 2 — n 
n 2 — 1 




1 1 1 

- + —y = — rr- 

n n z n + 1 
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6 A non- symmetrical qubit case for POVM 



In contrast to the previous section where we constructed an optimal POVM without 
any information on the state, in the upcoming part we will show some special, non- 
trivial cases, where some information is known. This is a much harder task, therefore 
the achieved results are not so strong, nevertheless, it is even more interesting than the 
fully unknown case. 

Let us assume, first of all, that we have a qubit: 6>i, 6 2 are not known, but #3 is 
given, so the density operator is only partially unknown. In this case the POVM will 
have three components {Ei, E 2 , E 3 }, the sufficient number of components to make the 
state estimation. Of course we can use a greater number of components but in that 
case the construction of the estimator is not trivial. 

We use the parametrization 

Ei = a (I + a ■ a), E 2 = b (I + b ■ a), E 3 = c (I + c ■ a), 

where a , b , c G K, and a, b, c are vectors in R 3 . From E\ + E 2 + E 3 = /, we have 
Oo + fy) + c o — 1 an d a a + b b + c c = 0. The positivity conditions are: 

< a ,b ,c \\a\\, \\b\\, ||c|| < 1. 



Similarly to the previous section we should minimize a determinant of the form A/B, 
where 



A — det V, B = det : 



ai a 2 
h b 2 



det 



(a*, a*) 
(a%b*) 



(a*,b*) 
(b*,b*) 



--: detC, 



V 



-pq 



-pq 
qbo 1 - q 2 



(a*, a*) (a*,b*) 
(a*,b*) (b*,b*) 



--.D-c-C 



p=l + a 3 9 3 , g = 1 + 636*3, a * = (ai,a 2 ), b* = (bi,b 2 ), 



where (3 is a constant depending on the particular unitarily invariant surface of the 
integration, < ft < 1, so c is a given constant on interval [0, 1 —9%]. 

For the minimizer we can assume some symmetry conditions: a = b = cq = |, and 
we will minimize A and maximize B independently to obtain the optimum: 

A = (du - c cu)(d 22 - c c 22 ) - (di2 - c Ci 2 ) 2 -> min, (16) 

B = C11C22 - (ci 2 ) 2 ->■ max. (17) 

Let us suppose first that a 3 and 6 3 are given (then the elements of D are constants) 
and we want to optimize the other variables. 

We know that 



(a + b, a + b) = cu + c 22 + 2c i2 + (a 3 + 6 3 ) 2 < 1. 



(18) 
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If C12 — 0, then from ffl8l) we have en + C22 < 1, hence 5 < 1/4. 

If C12 < 0, then B is maximal if Cyi is maximal, from fllSp we have an upper bound: 

1 - (03 + fr 3 ) 2 - en - c 22 
C12 < ^ • ( 19 ) 

Substituting this upper bound in ( IT7|) we have to maximize it in en and c 2 2- Using 
derivation we can obtain that it is maximal if en and C22 are maximal: 

c u = a \ + a\ < 1 - a 2 and c 22 = 6? + b\ < 1 - 6|. (20) 

Substituting this upper bound in ( FIT]) we get 

fl = !-ajj-a 3 &3-&|j. (21) 
which is optimal if a 3 = b 3 = 0. Then I? = 3/4, so it is a global optimum. 
Lemma 2 TTie following inequality is always true: 

di2 - c C12 < 

Proof. We have 

d 12 - c C12 = -(1 + 6303) (1 + 6" 3 & 3 ) - c(a x 6 a + a 2 6 2 )- 
Since c < 1 — #3 and 



axfei + a 2 &2 = (a*,b*) > -||a*||||6*|| > - af^/l - & 

so it is enough to show that 



(1 + 9 3 a 3 )(l + 9 3 b 3 ) > (1 - _ a 2^i _ 6 | 

The right-hand side does not depend on the signs, the left-hand side is minimal if 
@3 a 3 < 0, 9 3 b 3 < 0, so it suffices to prove for positive a 3 , 63, 6*3 that 



(1 - 9 3 a 3 )(l - 3 b 3 ) > (1 - ODsJl - a»yjl - bl 
This is true since from the Cauchy-Schwarz inequality follows that 



a 3 9 3 + a*y/l -9I<1 1 - a 3 9 3 > y/l - a*y/l - 91 

and a similar statement is true for b 3 . □ 
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Using this lemma, we get that A is minimal if C12 is maximal, and from there the 
solution is almost the same as in the previous case, we only have more complicated 
calculations. We substitute the upper bound ( TTTJT) in to (IT5|) and we can obtain that it 
is minimal if cn and C22 are maximal. Using the bounds ( 120]) we get for A a function 
of a 3 and 63; using differentiation we can obtain that 03 = 63 = gives the optimal 
solution here, too. 

So in both cases we have equality in (T2T)j) . so (a*, a*) = (b*,b*) = 1 and equality in 
PU), so (a\b*) = -1/2. Since ||a|| = ||6|| = 1, 

l El = l -{I + a-a) ) ^E 2 = ^(I + b-a) 

are projections. We have ||c|| = ||a + 6|| = 1 and this implies that 3E 3 /2 is a projection 
as well. 

Theorem 3 The optimal POVM for the unknown 61 and Q 2 can be described by pro- 
jections Pi (i = 1, 2, 3): 

3 

E * = l P *> ll Pi = \^ Tr ^: = 0, Tr P i P j = ~ for i^j. 



The optimal POVM is quasi-orthogonal to the subalgebra generated by 03 and sym- 
metric in the other directions. This result is in some sense the combination of Theorem 
1 and Theorem 2, with complementarity and symmetricalness both occurring in the 
same problem. The conditions of are valid, but with different constants than in the 
SIC-POVM case. 

It is an interesting question whether the existence of a conditional SIC-POVM in 
Theorem [3] is just a coincidence and it follows from the geometrical structure of M 2 (C), 
or there are such objects in higher dimensions. The following example solves this 
problem: 

Example 2 Let the elements of a POVM in M 3 (C) be 



E x = - 
7 



"1 


1 


r 




' 1 


£ 6 


e 2 ' 


> E3 = - 


' 1 


e 2 


e 3 ' 


1 


1 


1 


e 


1 


e 3 




1 


e 


1 


1 


1 




e 4 


1 


e 4 


£ 6 


1 



E d 



1 
7 



1 e 4 
: 3 1 



Ex — E-). Er — E*. En = Ea 



where e = exp (27ri/7). This POVM fulfills the conditions in ([6]) with constants k = 7, 
A = 7/3 and fi = 2/9. Moreover it is quasi-orthogonal to diagonal matrices, so it is a 
conditional SIC-POVM with respect to the diagonal matrices. □ 
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It is still an open question, however, that exactly in which cases a conditional SIC- 
POVM exists. Furthermore, does it minimize the determinant of the average covariance 
matrix? The answer to the first question is definitely a difficult problem, because it 
contains the question of the existence of SIC-POVMs. Our conjecture to the second 
question is that conditional SIC-POVMs are optimal if they exist, and so they can 
play an important role in quantum information theory, but an exhaustive investigation 
exceeds the scope of this paper. 

7 Discussion and conclusions 

A density matrix p e M n (C) has n 2 — 1 real parameters. If there are no known param- 
eters, several two-valued POVMs {F 1 ,I- F x }, {F 2 , 1 - F 2 }, ... , {F n 2_ u I - F n a_i} or 
a single POVM {Ei,E 2 , . . . ,E n 2} are considered. The estimator can be constructed 
as a linear function of the outcomes, and since they are multinomially and binomially 
distributed, their variance can be calculated easily. By averaging them on the invari- 
ant states and taking the determinant of the average covariant matrix, we will get 
comparable results. The optimum is at the quasi-orthogonal measurements (Theorem 
1) and at a symmetric informationally complete system (Theorem 2), respectively. If 
partial information is known, the concept is the same as previously, but as there are 
some known parameters, the domain of integration is restricted. From Theorem 1, it 
turns out that in the von Neumann case the optimal measurements are not only com- 
plementary to each other, but to the known subalgebras as well. For POVMs a qubit 
case is described in detail: conditionally symmetric informationally complete systems 
turn out to be optimal. The latter is a new concept, an example for the qutrit case is 
obtained too, further investigations are suggested. 
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